Super-Tonks-Girardeau gas of spin-1/2 interacting fermions 
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Fermi gases confined in tight one-dimensional waveguides form two-particle bound states of atoms 
in the presence of a strongly attractive interaction. Based on the exact solution of the one- 
dimensional spin-1/2 interacting Fermi gas, we demonstrate that a stable excited state with no 
pairing between attractive fermionic atoms can be realized by a sudden switch of interaction from 
the strongly repulsive regime to strongly attractive regime. Such a state is an exact fermionic ana- 
log of the experimentally observed super- Tonks-Girardeau state of bosonic Cesium atoms [Science 
325, 1224 (2009)] and should be possible to be observed by the experiment. The frequency of the 
lowest breathing mode of the fermionic super- Tonks-Girardeau gas is calculated as a function of the 
interaction strength, which could be used as a detectable signature for the experimental observation. 

PACS numbers: 67.85.-d, 03.75.Ss, 03.75.-b 



Introduction. — Exploring new quantum phases and 
understanding the striking consequences of correlation 
in strongly interacting atomic gases are at the frontier 
of current research in condensed matter physics and cold 
atom physics [l|, y|- In recent years, remarkable progress 
has been made in the experimental realization of funda- 
mental many-body model Hamiltonians, such as the Hub- 
bard model [Sl and the Tonks-Girardeau (TG) gas lJ,l5|i 
with unprecedented tunability. In general, attentions are 
devoted to the exotic properties of ground states and 
low-excited states. A recent experimental breakthrough 
is the realization of the super Tonks-Girardeau (STG) 
gas of bosonic Cesium atoms [6[ , which is a stable highly 
excited state of interacting Bose gas [7H9[ . In the exper- 
iment [6[, a one-dimensional (ID) Bose gas was initially 
prepared in the strongly repulsive TG regime, and then 
the STG gas was obtained by suddenly switching the in- 
teraction from strongly repulsive to attractive regime. A 
striking feature of the STG gas is its counterintuitive sta- 
bility against collapsing to its cluster ground state even 
in the presence of strongly attractive interactions. 



So far, the experimental study 
oretical works 



and most of the the- 



12l | on the STG gas have focused on 
the bosonic system. In this work, we study the possible 
realization of the Fermi super Tonks-Girardeau (FSTG) 
state in a ID Fermi gas. As the ID Fermi gas with tun- 
able interaction strengths has already been experimen- 
tally realized 13| , it is promising to directly observe the 
STG state in ID Fermi gases. Stimulated by the exper- 
iment of the Bose STG gas, we first prepare a strongly 
repulsive spin-1/2 Fermi gas, and then suddenly switch 
the interaction to the strongly attractive regime. By this 
way, we can access a stable highly excited state of the at- 
tractive Fermi gas which does not fall into its attractive 
ground state. In the strongly attractive limit, atoms with 
different spins form tightly bound fermion pairs |14|-|l6|. 
It has been shown that the ground state of a ID attrac- 
tive spin-balanced Fermi gas is effectively described by 
the STG state of bosonic pairs of fermions with attrac- 



tive pair-pair interaction [16[. The FSTG state being 
studied in the present work is essentially the lowest gas- 
like excited state composed of unpaired fermions which 
is totally different from the bosonic STG state composed 
of tightly bound fermion pair proposed in Ref. [l6| . 

Interacting Fermi model. — We consider a system of 
N = iV-|- + N^ spin-1/2 fermions in a tightly confined 
waveguide described by the effective ID Hamiltonian 
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where gij. — — 2?i /{maid) ~ ^ c/m is the effective ID 

interaction strength and aio the effective ID scattering 

length [17[. Without loss of generality, we assume that 

N\, < N^. The ID interacting spin-1/2 Fermi gas is only 

solvable for the homogenous case [IJ]. However, in the 

infinitely repulsive limit, a generalization of Bose-Fermi 

r— T I — I I — i 

mapping [18| to the spin-1/2 Fermi system [19|,l20[ allows 

us to construct analytically exact solution of ID Fermi 
gases even in trap potentials. In this work, we focus on 
the homogenous system which can be exactly solved by 
the Bethe-ansatz (BA) method. The trapped system can 
be studied in the scheme of the local density approxima- 
tion (LDA). 

The model ^ is exactly solved by the BA method [1^ 
with the BA wavcfunction 

tp{xi ,■■■ ,xn) = E E ^(^Qi < ■■■ < xqn) X 
Q P 

N 

[Q,P]exp[i^fcpjXQj], (2) 

where ki represent quasimomenta, P and Q represent 
permutations of ki and xi, respectively. For the eigen- 
state with the total spin S = N/2 - M {M = A^j,), the 
coefficient [Q, P] can be explicitly expressed as [Q,P] = 

Et=i '^ivTi ,yT2,- ■ -vtm ; P) ]\j=i XvT^ 4 YIx^^vt- X:r.,t, 



where Xxi,\ [Xvja) denotes the up (down)-spin, T is 
a combination of M down-spins in N particles, {ht } 
are M elements of T{xi\, and '^{yTnUT-i^ ' ' 'HTm^P) = 
Y^j^ A{R) rij-i Fp{^Rj 1 TJTj ) with R being the permuta- 
tions of As, A(R) = e(i?)n„^^;(AR, - A^^, - zc), and 



A{R) = e(i?)n,<i(Afl. 
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Aji- — ic/2). The parameters kj and A „ ar e determined 
by the Bethe-ansatz equations (BAEs) [l 
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y 2tan-i(l^^i— ^) = 2^J„ + 2ytan-i( ^""^^ ). 
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Here 



eigenenergies are given by £■ = ^ ^^. ,„^ . 
both kj and A^ are real numbers if c > 0. The ground 
state solution corresponds to Ij = (iV+l)/2— j and Ja = 
(Af-f l)/2— a. In the limit of cL ^ 1, A^ arc proportional 
to c, but kj remain finite, therefore the quasi-momenta 
can be given approximately 



k,L = 27rI,-A + Oi\c\-'') 



(5) 



with ( ~ X]a=i /A / yi+i/A - It follows that the ground 
energy in the strongly repulsive limit reads 
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which is consistent with the result in Ref . [21| up to order 
of c~^. In the limit of c — > oo, the ground energy is 
identical to that of a polarized N-fcrmion system. 

FSTG state. — If the interaction is attractive, the 
ground state is composed of iV — 2M real ki and 2Af 
complex ones. In the limit —cL ^ 1, the complex 
solutions take the 2-string form: fc^ ~ A^ + |i, and 
kM+a ~ Aq, — |i. Except the complex solutions, the 
BAEs also have real solutions for c < 0, which, however, 
correspond to some highly excited states of attractive 
Fermi systems. The FSTG state corresponds to the low- 
est real solutions of BAEs ^ and (|3]) with c < 0. In this 
case, Aq, go infinite and kj remain finite with \c\L — > oo, 
thus the momenta are given by 



k,L = 2nI,+Cj^ + 0{\cr'). 



(7) 



Despite the C in Eq. ([7]) having the same form as in 
Eq. (O, generally C,{c) ^ C(~c) since the solutions Aq, 
of Eq. (J4|) are not symmetric for c and — c. However, in 
the strong coupling limit, up to order of c~^ Eq.Q be- 
comes 2Ntan-\^) = 2-KJa + 2 Y!p=i tan-i(^^^i^) + 
0(|c|~^), which is invariant under the operation P : {c -^ 
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FIG. 1: (color online) (a) Quasi-momentum distributions for 
tlie ground state of the repulsive Fermi gas and the FSTG 
state of the attractive Fermi gas with different values of 7. 
(b) The energies Eftg (stars) and Efstg (dots) vs 7. 



— c, Aq, — > — Aq}. Therefore we have C,[c) = ({—c) up to 
the order of c~^. The energy of the FSTG gas in the 
strongly attractive limit is thus given by 
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In the limit of \c\ -^ cxd, we have Epsra — EpTo and 
kj = /j27r/L for both the Fermi TG and the FSTG gas. 
In Fig. 1(a), for an example system with A'^ = 10 and 
M = 5, we show the BAE solutions of kj for the repul- 
sive Fermi gas and the attractive FSTG gas with differ- 
ent values of 7 = c/p, where p = N/L is the particle 
density. The quasimomentum distributions for the re- 
pulsive Fermi gas and the STG gas approach the same 
limit from different sides when I7I goes infinite. Corre- 
spondingly, EpsTG and EppG Silso approach the same 
limit Einf = ^j^N{N'^ - 1) as shown in Fig.l(b). 

The FSTG state can be achieved through a similar 
sudden switch as in Ref. [6| . The system is first prepared 
at the ground state in the strongly repulsive regime, i.e., 
\'^{t = 0)) = |(/3o(c)). After a sudden switch into the 
strongly attractive regime with c' < 0, the wavefunction 



is given by |^(t)) 



-iHt 



Vo{c)) = Y^-e 



iEit 
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where '^i{c') is the i-th eigenstate of the Hamiltonian 
with parameter c' and ai ~ {Lpi(c')\fo{c)) . The proba- 
bility for the system staying in a state \ipi{c')) is given 
by |aip = |((/5i(c')|(/'o(c))p. We note that the wave- 
functions ifsTGic') and 'Po(c) are identical when c' = —00 
and c = 00, and thus one can expect the probability 
of the system transforming from the Fermi TG gas to 
STG phase to be close to 1 for large \c'\ and |c|. In 
Fig 2, we display transition probabilities from the ini- 
tial ground state with c > to the STG phase with 
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FIG. 2: Transition probabilities from the Fermi TG gas to 
FSTG piiase for systems with N ^3,A, 5. 



FIG. 4: Energy vs 7 for states with different total spin S. 




FIG. 3: (color online) Local correlation energy vs 7 for the 
FSTG state and ground state of the repulsive Fermi gas. 



c' ~ — c for difFcrcnt-size systems. Our results show that 
the overlap of wave-functions approaches the limit of 1 
for large I7I, although it decreases as the system size 
increases. In the strongly interacting limit, the transi- 
tion probabilities can be approximately represented as 
P{N,M,-f) = 1 - a{N,M)/j^. The upper bound of 
the parameter a{N, M) can be estimated as a{N, M) < 
647r^7VAf ^ through the expansion of wavefucntions to or- 
der of 1/c. The transition probability for a larger sys- 
tem is expected to approach 1 if I7I ^ MN^/"^. How- 
ever, the calculation for a large system becomes a very 
time-consuming task due to the calculation of multidi- 
mensional integrals. 

In thermodynamic limit where iV, M, and L go infinite 
but p = N /L and m = M/L remain finite, the BAEs can 
be expressed in the form of the coupled integral equations 



IJ). The energy per particle reads 5(7) = ^P^e(7). For 



the case with m = p/2 = J_ cr(A)(iA, C, can be calcu- 
lated via the integral C, = /^ dxAa{x) / [x^ + 1) = 2 In 2, 
where ct{x) = l/[4cosh(7ra;/2)]. Therefore we can ob- 
tain the energy expansion 657^0(7) = 'xl-'^ 4-41n2|7|~^ + 



12(ln2)2|7|-2] + 0(7-3) for large |7|. The energy in 
the whole area can be calculated by numerical solving 
the coupled integral equations. Comparing the expan- 
sion result with the exact numerical result, we find that 
they agree very well in the regime of large 7, for ex- 
ample, \C^xpansion ^nu7n\/ \^num\ ^ J^U lOr 7 — ^^■ 

A characteristic of the STG gas is that it has stronger 
correlations than the TG gas. In the regime of I7I 3> 
1, the local two-particle correlation function (72(7) = 
p^de{'y)/d'y can be directly derived from 6(7), which 
gives g2il)TG/p^ ~ (47rV3)(ln2|7|-2 - 6(ln2)2|7|-3), 
and ff2(7)sTG/p' ~ (4^V3)(ln2|7r2+6(ln2)2|7|-3). 
Thus we have 52 (—|7 1 )sTG > 52(|7|)tg asshowninFig.3. 
In terms of terminologies of Tomonaga-Luttinger liq- 
uid (TLL) theory [22], the strongly repulsive phase of 
the spin-balanced Fermi gas corresponds to a TLL with 
the charge TLL parameter K^ « (1 -h 41n2/|7|)/2 > 1/2 
22l |. The FSTG phase corresponds to a highly excited 
gas-like state where unpaired particles are strongly cor- 
related. This strongly collective behavior may be phe- 
nomenologically described by Kc ^ {1 — 41n2/|7|)/2 in 
the strongly interacting limit, which is smaller than 1/2. 
For both cases, the spin TLL parameter K^r = 1 due to 
spin-rotational invariance. 

Degeneracy of the FSTG state. — In comparison to the 
Bose system, the ground state of a spin-1/2 system is 
highly degenerate in the TG limit due to the fact that 
states with different total spins have the same energy 
|20[. However, for a large but finite interaction strength 
the degeneracy is broken and the true ground state is 
the state with the lowest 5*. For the spin-1/2 system 
described by ([T]), one can understand this fact from the 
energy expression ([5]). The term ( is M-dependent (M = 
N/2-S) and we have CiMi) < ((^^2) if Mi < M2, which 
leads to Eftg{S2) < EFTciSi) for 5*2 < ^i. This is 
consistent with the Lieb-Mattis theorem [23| . The energy 
difference is proportional to 1/c and vanishes as c — )■ 
00. On the other hand, for the FSTG state we have 
Efstg{S2) > Efstg{Si) for ^2 < Si according to Eq. 



([8]), i.e., the FSTG state with the smaller S has higher 
energy. To give an example, we calculate the ground state 
energy and the energy of the FSTG state for a system 
with A^-|- = N^ = 2. As shown in FigUJ energies for states 
with different total spins approach the same limit of the 
polarized Fermi gas as \c\ — > c«. 

Experimental detection. — To realize the FSTG gas, 
one can first tune the interaction of the ID Fermi gas to 
the strongly repulsive regime by the Feshbach resonance 
13| , and then suddenly switch the interaction across the 
resonance point. Similar to the bosonic case, one can 
measure the frequency of the breath mode of the FSTG 
gas subjected to a weak harmonic confinement along the 
axial direction, which is sensitive to various regimes of 
interaction. For the Fermi gas in a harmonic trap with 
Vext = muj'^x^ /2, we can determine the density distribu- 
tion of the STG gas within the LDA. According to the 
LDA, the system is in local equilibrium at each point 
X in the external trap. The density distribution of the 
FSTG gas is then obtained via the local equation of state 
Mo = l^[p{x)] + Vext{x) [ii-lii]. Here /^(p) = dp[pe{p)] is 
the local chemical potential with e(p) the energy den- 
sity of the homogenous FSTG gas to be determined by 
numerically solving integral BA equations, and yUo is de- 
termined by the normalization condition / dxp(x) = N . 
Following Rcfs. [25|, |26|, we calculate the frequency of 
the lowest breathing mode from the mean square radius 
of the trapped FSTG gas via a;^ = -2{x'^) /{d{x'^) / dul) 
with {x^) = J p{x)x^dx/N. The solid line in Fig. 5 shows 
the frequency of breathing mode of the attractive FSTG 
gas as a function of the interaction strength Na\^/a^ 
with the harmonic oscillator length ax = yh/rnujx. The 
frequency of the breath mode for the FSTG gas exhibits 
a peak with a maximum of w^/w^ about 4.3. We also give 
results of the repulsive Fermi TG gas (the dashed line) 
and the ground state of the attractive Fermi gas (the 
dotted line) for comparison. These results are essentially 
based on the sum-rule approach and provide generally an 
upper bound on the frequencies [25[. 

Summary. — In summary, we study the realization of 
the FSTG phase in the interacting spin-1/2 Fermi gas. 
Starting from the ground state of a strongly repulsive 
Fermi gas, the FSTG state can be realized by a sudden 
switch of interaction to the strongly attractive regime. It 
is shown that the FSTG state is stable against forming 
pairing states even in the presence of the strongly attrac- 
tive interaction between fermionic atoms with opposite 
spins. We also calculate the lowest breathing mode fre- 
quency of the FSTG gas which may be detected by the 
experiment. 

Note added. — While this work was being prepared for 
submission, a related manuscript appeared [27| , in which 
the FSTG gas in a harmonic trap is studied. 
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FIG. 5: Square of the lowest breathing mode frequency vs the 
interaction strength Na\^/a^ for the FSTG gas (solid line), 
the repulsive Fermi TG gas (dashed line) and the ground state 
of the attractive Fermi gas (dotted line). 
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